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ABSTRACT 

The hypothesis that massive stars form by accretion can be investigated by 
simple analytical calculations that describe the effect that the formation of a 
massive star has on its own accretion flow. Within a simple accretion model that 
includes angular momentum, that of gas flow on ballistic trajectories around a 
star, the increasing ionization of a massive star growing by accretion produces 
a three-stage evolutionary sequence. The ionization first forms a small quasi- 
spherical HII region gravitationally trapped within the accretion flow. At this 
stage the flow of ionized gas is entirely inward. As the ionization increases, the 
HII region transitions to a bipolar morphology in which the inflow is replaced 
by outflow within a narrow range of angle with about the bipolar axis. At 
higher rates of ionization, the opening angle of the outflow region progressively 
increases. Eventually, in the third stage, the accretion is confined to a thin 
region about an equatorial disk. Throughout this early evolution, the HII region 
is of hypercompact to ultracompact size depending on the mass of the enclosed 
star or stars. These small HII regions whose dynamics are dominated by stellar 
gravitation and accretion are different than compact and larger HII regions whose 
dynamics are dominated by the thermal pressure of the ionized gas. 

Subject headings: stars: early type — stars: formation — ISM: HII regions 



Introduction 



The hypothesis that massive stars form by accretion can be investigated by simple 
analytical calculations that describe the e ffect that th e formation of a massive star has on 
its own accretion flow. An earlier paper (IKetol l2002bf) stud ied the effect of star formation 
in a spherically-symmetric, steady-state flow (IBondil Il952l ) and showed how a molecular 
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accretion flow may pass through an ionization front and continue toward the star as an 
ionized accretion flow. In this case, the HII region is the ionized inner zone of a continuous 
two-phase accretion flow. More generally, angular momentum in an accretion flow will result 
in a flattened geometry with different consequences for the structure of the ionized accretion 
flow and HII region. 

This paper investigates the effects of the ionizing radiation of a massive star on an 
accreti on flow with angular momentum, modeled as gas flow on ballistic trajectories around 
a star (jUlrichlll976l ). The model is described by a few parameters: the gas density in the flow, 
the angular momentum, the mass of the star, and the flux of ionizing radiation. Alternatively, 
the model can described in terms of three characteristic radii: characteristic radii, Ri, the 
radius of ionization equilibrium, R b , the radius where the escape velocity equals the sound 
speed, and Rd, the radius of disk formation where the gravitational and centrifugal forces 
balance, 

Models with different ratios of these three radii produce HII regions with different mor- 
phologies and dynamics. In particular, increasing the ratio Ri/ Rb radius of ion ization results 
in a progression of HII regions from quasi-spherical, gravitationally-trapped iKetol (j2002bl . 
20031 ) through b ipolar to a final morphol o gy simi l ar to t hat around a photo-ev a porating disk 
as des c ribed by iHollenbach et al.l (1994 ) . lYorkel (11995). lYorke fe Weld (1996). Ilizano et al. 



(Il996h . I Johnstone. Hollenbach. fe Ballvl (119981 ). and lLugo. Lizano fe Garavl (120041 ) . This pro- 
gression suggests an evolutionary sequence for an HII region around a star that is growing by 
accretion. As the star gains mass, the location of the radius, Rb, will increase proportionally 
with the mass of the star. The radius of ionization equilibrium, Ri, will also increase but as a 
higher power of the mass because the flux of ionizing radiation depends on the temperature 
of the star. Thus this simple model of ionization within an accretion flow provides an under- 
standing of the relationship between the several different observed HII region morphologies 
on the hypercompact and ultracompact scales. 



2. Relevant results from the spherical model: the 3 evolutionary stages of an 

HII region 



Stellar structure models suggest that massive stars that form by contraction have a 
shorter pre-main sequence phase (PMS) than low mass stars, and that massive stars that form 
by accretion h ave no PMS phase at all because massive stars begin hydrogen burning while 
still accreting JPalla fe Stahler!ll993 : Il3eech fe Mit alasll 1 9 94j: I Chieffi . Staniero fe SalariJ ~ 



1995 



Bernasconi fe Maeder 111996: Behrend fe Maeder 



200ll : iNorberg fe Maeder II200CH : iKeto 



2003h . 



These calculations suggest that a massive star that is forming by accretion will reach the 
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main sequence with a mass below that which would produce sufficient ionizing radiation 
to maintain an HII region around the sta r. The non-existent HII region at this stage may 
be described as quenched (jWalmsleyl Il995l ) . Because the production of ionizing photons in- 
creases as a higher power of the stellar mass than does the accretion rate, eventually the flux 
of ionizing radiation, J*, fro m a star th at is gaining mass by accretion, will exceed the flux 
of neutral gas onto the star (jKetol 120031 ) . 



T. 



where r* is the stellar radius, n# is the the number density of neutral gas, and v the inflow 
velocity. 

If the density of the accretion flow decreases with distance no fa ster than r~ 3 / 2 , then the 

HII r egion that develops in the flow will be bounded at a distance, Ri (IFranco. Tenorio-Tagle fc Bodenheime 
19901 ). defined by the balance of the ionizations and recombinations (ISpitzerlll978l . eqn. 5.20), 



J $ 



Airr n e {r)a dr = 



(2) 



This equation incorporates the "on-the-spot" approximation with a the recombination rate 
to all Rydberg levels of H above the first level, n — 1. 

The small HII region that develops just as the quenched phase (equation [1]) ends cannot 
expand hydrodynamically because the outward pressure of the ionized gas is less than the 
inward gravitational force of the star. Instead, the ionized gas forms part of a continuous 
accretion flow with an outer molecular phase and an inner ionized phase separated by a 
static ionization front within the flow. The HII region, which at thi s stage is the inner part 
of the accretion flow, may be described as gravitationally trapped ( jKetol l2002bl . 120031 ) . As 
the star continues to gain mass and its flux of ionizing radiation increases, the boundary of 
the HII region, whose location is determined by ionization equilibrium rather than pressure 
balance, will be found at greater radii. Once a radius is reached where the escape velocity 
from the star is below the sound speed of the ionized gas, the HII region will begin to expand 
hydrodynamically. This radius is approximately equal to the Bondi-Parker trans-sonic radius 
where the velocity of the incoming accretion flow first reaches the sound speed of the ionized 



gas, 



R b = GM*/2c 2 



(3) 



where M* is the mass of the star and c is the sound speed. The evolution of the HII region 
then transitions to a phase of pressure- driven expansion in which the gravitational attraction 
of the star is negligible over most of the HII region. These pr essure dominated HII regions are 
adequately described by models with no gravitational force ( )Spitzerlll978l ; iDyson fc Williams 
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19801 ; IShulll992l ). The model for the development of an HII region within a spherical accretion 
flow thus suggests three stages: 1) quenched or non-existent, 2) gravitationally trapped, and 
3) pressure driven-expansion. The divisions between the stages occur when the radius of 
ionization equilibrium, R iy (equation [2]). 



3. A non-spherical accretion flow 



One simple model of accretion that includes rotation equates the streamlines of an ac- 
cretion flow with ball i stic trajectories around a poin t mass ( )Ulrichlll976l ; ICassen fc Moosman 
198ll ; IChevalierlll983l ; iTerebey. Shu fc Cassenlll984l ) The flow is determined by the gravita- 
tional attraction of the point source subject to conservation of angular momentum and mass. 
In particular, the model ignores the pressure and self-gravity of the gas. The initial distribu- 
tion of angular momentum is that of solid body rotation on an arbitrary radius or conceptual 
spherical surface some distance from the star. Thus on each trajectory, the specific angular 
momentum, T oc sin# ? where 9 is the initial polar angle. The gas density is described 
by the conservation of mass along streamlines. In this model, an accretion disk develops 
as the gas density increases along inwardly converging streamlines flattened by the angular 
momentum of the flow. A limitation of the model is that the mid-plane density is not defined 
whereas a more com plete model would describe the density across the mid-p l ane as a pres- 



sure supported disk (IShakura fc Sunvaevl Il973l ; iLynden-Bell &: Pringld Il974j ; iPringld Il981 



HartmanrJfl9~9~8l : Iwhitnev et al.ll2003h 7 This deficiency may be ignored if, as in the examples 
that follow, the scale of the HII region developed within this accretion flow is larger than 
the scale height of the disk (see S|6]). The equations for the velocity and density at a position 
(r,9) are Juirichlll976h . 

*=-{—) ' 



COS 9 \V2 



COS I 



V<f, 



COS #0 

GM\ 1 / 2 sin 9 
sin 9 



1 + 



COs6> 

COS 9q + COS #o 



cos 6*o sin 9 
cos9\ 1 / 2 



1/2 



cos 9r 



(4) 
(5) 
(6) 



where, 



Rd cos9 sin 9 2 



(7) 



cos #o — cos 9 

This last equation includes the parameter, R d , which is the radius at which the gravitational 
force equals the centrifugal force in the mid-plane, 



GM/rl 
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This radius, approximately where the accretion flow transitions from a quasi-spherical inflow 
to a rotationally supported disk, may also be written in a form analogous to the Bondi-Parker 
transonic radius as, 

R d = GM/v 2 k (9) 

where t^, the orbital veloci ty at Rd, replaces the sound s peed in equation [31 The gas density 
from mass conservation is (IMendoza. Canto. Ragall2004l ). 



n 



UqT 



-3/2 



1 + 



cos0\-V2 



COs6> 



l + r-^cos 2 ^-!)) -1 



(10) 



where the number density, n , is defined through the mass density, fin Q , by the mass accretion 
rate, 

M = fin A7iR 2 d v k (11) 



4. The development of an HII region in the non-spherical accretion flow 

An HII region will develop at the center of the flow defined by equations H] through [IT] 
once the ionizing flux of the accreting star satisfies the condition in equation [TJ However, 
in a non-spherical flow, the gas density is not radially symmetric, and therefore the HII 
boundary, Ri, is at larger distances in directions where the gas density is lower. In the 
axially symmetric model above, the boundary is therefore a function of the polar angle, 9. 
The spherical model of $2] suggests that if Ri > Rb the HII region expands hydrodynamically. 
Thus while somewhat idealized, in the non-spherical model the HII region may be expanding 
over a range of angle 9 about the polar axis while the accretion flow continues to flow into 
and through the HII in a different range of angle 9 around the plane of the disk. 

This simple model of the ionization of an accretion flow does not explicitly include the 
effect of a stellar wind on the HII region. We know that unembedded early-type stars produce 
powerful winds that at their terminal velocities have mechanical energies, pvf, t hat exceed 



the mechanical energies of star-forming accretion flows (lLamers fc Cassinellilll999l ) . However, 
because the winds accelerate off the surface of the star the wind energy is negligible near 
the star and, in the idealized spherical model, the wind may be suppressed by the accretion 
flow. Once inflow ends because the gas within the HII region begins expanding (Ri > Rb), 
the wind may be expected to break out and dominate the hydrodynamics. Thus in the 
non-spherical model, the following approximation is adopted. At a polar angle, 9, where 
Ri(9) > Rb, the inflow is replaced by the simplest model of a stellar wind (Parker 1958). 
This wind solution is the analogue of the Bondi accretion flow and described by the same 
equation, 

d fv 2 \ dP GM 



*ItJ+*+ — = ° (12) 
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but with the boundary conditions that the flow is subsonic at the stellar surface and su- 
personic beyond the transonic radius. In the case of non-spherical flows, the location of 
the transonic radius depends on the geometry of the fl ow, specifically on the divergence of 
the streamlines about the star (IKopp Holzerl Il976l ). However, if the divergence is ap- 
proximately spherical, the transonic point will be approximately located by equation [31 For 
simplicity, this approximation is adopted. The resulting model of a wind in the polar di- 
rections and accretion at lower latitudes is inconsistent in that the wind solution is driven 
by hydrodynamic pressure, but pressure is neglected in the accretion flow. Nonetheless, 
the examples in §6] will show that this approximation is useful on a conceptual level and 
produces models that describe the morphology of observed HII regions. Different model s 
that describe the interaction of a w ind and accretion flow inc lude IWilkin Stahlerl (119981 ) , 



Mendoza. Canto. Ragal (120041 ) , and ICunningham et al.l (120051 ) . 



5. The 3 R's of hypercompact HII regions 

The equations in sections £|2] and 2] fully describe the model of an HII region in a 
non-spherical accretion flow. The morphology of the HII region depends on a few physical 
quantities that are specified as parameters of the model: the gas density and angular momen- 
tum in the accretion flow, and the mass and flux of ionizing radiation of the star. Because 
these physical quantities have different units, it is helpful to organize them into parameters 
of the same physical dimension as 3 characteristic radii. The radius of disk formation, Rd 
(equation EJ), along with the Bondi-Parker radius, Rb (equation [3]) , and the radius of ion- 
ization equilibrium, Ri (equation [2]), may be called the 3 R's of HII regions in that they 
describe the basic structure of the accretion flow and HII region. 

We may also combine these 3 characteristic radii into non-dimensional ratios, for ex- 
ample Ri/Rd, Ri/Rb, Rd/Rb to see that the models are scale-free, and that it is the relative 
magnitudes of the 3 radii that are important in describing the structure of the accretion flows 
and HII regions. Thus a model appropriate for a common accretio n flow onto several O stars 



such as in G10.6-0.4 with a total stellar mass of a few hundred M ( jKeto fc Woodll2006l) may 



have the sa me structure as a model for a single B star such as IRAS20126 (ICesaroni et al. 



19971 . Il999i ) provided that the three ratios are the same. 



Because the model structures depend on only 3 characteristic radii or 3 ratios, then only 
a few model examples are required to illustrate the possible variations as well as suggest an 
evolutionary sequence for HII regions. So that these models may be more easily related to 
the densities and stellar types familiar from observations of high mass star forming regions, 
the examples in the next section are presented in units appropriate for early B stars and late 
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O stars since these are most commonly observed massive stars. These models are applicable, 
i f scaled up, to accret ion flows as large as those onto groups of early O stars such as G10.6-0.4 
fceto fc Woodll2006h . 



6. Examples and evolution 
6.1. Ionization 

The first example illustrates the effect of increasing the ionizing flux within the same 
accretion flow. This comparison illustrates the differences as Ri/Rb is greater than, approx- 
imately equal to, or less than one. (In this example for the purpose of this comparison, we 
set Rd ~ Rb.) As the ionizing radiation from the star increases, and Ri increases relative 
to Rb, the HII region transitions through three phases from nearly spherical inflow of the 
ionized gas through bipolar to nearly spherical outflow. 

The initial model is comp osed of a Bl star with a mass o f 18 M and a flux of ionizing 



radiation of 10 45 photons s 1 (IVacca. Garmany fc Shulllll996l ) within an accretion flow with 



R d = 42 AU (r = 0.04 kms" 1 pc), R b = 43 AU, and a gas density, n of 10 7 cm" 3 . The 
mass accretion rate (equation [11]) is 6 x 10~ 6 M yr _1 . Stellar structure calculations suggest 
that this ac c retion rate is near th e minimum required for the formation of stars of this mass 
ket J l2003l : iKeto fc Woodl 120061 ). Higher rates might be required to form stars of higher 



mass. 

Figured] (top) shows that at this flux level and density, although the HII region has a 
slightly bipolar shape owing to the flattening of the flow, the boundary of the HII region 
is entirely with the Bondi-Parker critical radius (Ri < Rb at all polar angles). Thus the 
HII region is trapped within the gravitational field of the star, and accretion may proceed 
through the HII region at all angles. 

Increasing the stellar mass to 20 M Q and the flux of ionizing radiation to 7 x 10 45 photons 
s -1 , corresponding to an BO. 5 star, increases the extent of the ionization sufficiently so that 
Ri > -Rb in a narrow range of angle about the bipolar axis. Within this angular range, we 
assume that the HII region begins to expand and the structure and dynamics are described 
by an outflow or wind driven by the thermal pressure of the ionized gas. Accretion through 
the HII region continues at lower latitudes (figure [I] middle) . 

Increasing the stellar mass to 22 M Q and the flux of ionizing radiation to 10 47 photons 
s -1 , corresponding to an 09 star, results in a model with Ri > R b , and thus outflow, 
everywhere except at a narrow range of angle just around the disk (figured] bottom). This 
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third- stage si 


;ructure is similar to the photo-evaporating disk model ( 


Hollenbach et al. 


1994 


Yorkc 


1995 




Yorke & Welz 


1996: 


Lizano et al. 


1996; 


Johnstone. Hollenbach. & Ballv 


1998 



Lugo. Lizano &: Garayll2004l ) except that here the outflow is radially from the center of the 
HII region rather than vertically off a disk. A more complete model would include the 
evaporation off the disk in the outflow. 

In this scenario, the bipolar phase is relatively brief in the evolution between the gravi- 
tationally trapped and outflow phases. In this example, the bipolar phase corresponds to a 
range of stellar mass between 20 and 22 M Q . 



6.2. Angular momentum 

The second example shows the effect of different values of the angular momentum on the 
structure of the accretion flow. This example illustrates the difference as Rb/Rd is greater 
or less than one. Flows that have relatively low angular momentum will produce what may 
be described as a fat accretion torus around a growing star while those with relatively high 
angular momentum will produce a thin disk. Because the density of the flow increases inward 
following mass conservation regardless of the angular momentum, the difference between the 
two cases relates to where in the flow the gas density becomes observationally significant 
with respect to the flattening of the flow. The flow becomes observable as the gas density 
approaches the critical density for collisional de-excitation, n c = Aij/Cij, where Ay is the 
Einstein A coefficient (sec -1 ) and CV,- is the collision rate (cm -3 sec -1 ) for the particular 
molecular tracer employed in an observation. In a model with low angular momentum, the 
gas will reach the critical density before flattening into a disk, and the molecular accretion 
flow will be detected as a quasi-spherical rotating flow or fat torus. In the high angular 
momentum case, the flow will form a disk before the gas density in the surrounding flow 
exceeds the critical density. In this case the observations will tend to see the flattened disk. 

Figure [2] shows two variations of the same model. Both models include a star of 20 
M and flux of ionizing radiation of 3 x 10 46 photons s -1 corresponding to type BO - 09.5. 
At this mass, Rb = 54 AU. In one model, the angular momentum is Tq = 0.04 kms -1 
pc, while in the second model the rotation rate is 2 times higher. The values of Rd are 
therefore 38 AU and 152 AU respectively. Because of the greater inward increase in the gas 
density in the flattened, high-angular-momentum flow, the density in this flow is set lower, 
uq = 1 x 10 6 cm -3 , rather than uq = 3 x 10 7 cm -3 , so that similar densities are maintained 
at the center of the flow, and therefore, similar opening angles for the bipolar outflow. The 
mass accretion rates for each model are 6.6 x 10 -6 and 1.8 x 10 -5 M yr -1 . The two density 
distributions shown in figure 2 suggest that accretion disks and "fat toroids" may be different 
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expressions of the same model accretion flow with different values of the angular momentum 
or equivalently different values of Rd with respect to R^. 



6.3. Limitations of the model 



The model of accretion on ballistic trajectories in section §3] does not include pressure 
forces and therefore does not describe a pressure supported disk at the mid-plane. The 
density in such a disk could potentially affect the morphology of an HII region depending 
on the structure of the disk. Little is known either observationally or theoretically about 
disks around massive stars. If these disks are similar to those around low mass stars, then 
we can calculate how a scaled-up, low-mass disk would affect the HII region morphologies in 
the examples above. To do this we add to the model a ccretion flows in the examples above 
the density of a disk described as (jWhitney et al.ll2003l equation 3), 



p[r) 




exp 



2 \H r, 



(13) 



where r is the radius in the disk mid-plane, z is the height above the plane, a = 2.25, 
and the scale height, H(r) = H (r/R*) p with H = 0.1i2* (AU) and (3 = 1.25, The density, 
Po is defined by the assumption that the total mass of the disk within Rj is 0.1M*. These 
additional calculations (not shown) demonstrate that the inclusion of such a disk does not 
affect the morphologies of the HII regions in these examples. This is because the scale height 
of the pressure-supported disk at the radii where the boundary of the HII region meets the 
mid-plane is small compared to the thickness of the "disk" created by the flattening of the 
accretion flow. Because this pressure-supported disk is an arbitrary addition to the model 
that does not affect the results sought in this investigation, this disk is not considered further. 



7. The importance of stellar gravity and accretion on Hypercompacts and 

small Ultracompacts 



The morphologies produced by models of HII regions that develop within accretion flows 
are potentially quite varied and match many of the observed morphologies, particularly those 
observed at the smallest scales. In particular, the bipolar HII regions are predominantly a 
sm all scale phenomenon . For e xam ple, there is no cla ssification for bipolar morphology in 



the Wood fc Churchwe 



ll (119 89) and iKurtz et al.1 (119941 ) survey of ultracompact HII regions 



whereas iDepree et al.l (120051 ) find it necessary to introduce this class to describe the mor- 



phologies of hypercompact HII regions. 
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Th e size scale, ~ 0.01 pc (2000 AU), of a large hypercompact or small ultracompact HII 
region (jKurtzl l2000l ) . is about 20 times the radius, Rb around a single massive star. What 
then is the relevance of the stellar gravity for most observed HC and UC HII regions? 

First, O stars do not appear to form alone, but in small groups or clusters of early type 
stars that may also contain some number of lower mass stars. The radius, Rb, scales with 
the attracting mass. Groups or small clusters of stars with a few hundred M contained 
within a common HII region are suggested by observations of a few o f the brightest HII 
regions that show molecular or ionized accretion on scales of 10 s AU (IZhang 



k Ho 1997 



Young. Keto fc Hoi Il998l : iKetcl l2002al : ISollins et al. 1 120051 : iKeto fc Wood! 120061 ). Thus the 
gravity of groups or small clusters may be significant even on the size scales of ultracompact 
HII regions. 

Second, the HII region need not be strictly smaller than Rb for the stellar gravitational 
attraction to affect its structure. For example, if the HII region is in effect an expanding 
wind, as in the model s hown in figured] (bottom), then because the wind solution is quasi- 
hydrostatic (jKetol 120031 . equation 4) inside of Rb, the HII region will be quite dense within 
this radius even though the boundary, > R b . Outside of Rb, mass conservation within 
the converging accretion flow and within the diverging outflow generally requires that the 
flows have density gradients. Thus even beyond Rb, in the region dominated by pressure 
forces, the density is not uniform. High frequency observations would show a bright core 
whereas low frequency observations of the same HII region would show surrounding low-level 
extended emission. 



8. Jets, outflows and disks 



The outflows described here are simple models of ste llar winds driven by the pressure 
of the ionized gas (Parker isothermal wind (jParkerl Il958l )) in the HII region and confined 
in a specific way defined by the ratio of Ri/Rb as a function of angle, but nonetheless 
essentially by the geometry of the surrounding accretion flow. Thus these outflows require 
both significant ionizing flux from the star and a massive accretion flow. In this respect they 
are quite different from the more familiar bipolar jets or outflows associated with young low- 
mass stars. These latter are not completely understood, but are thought to be driven by the 
twisting of magnetic fields by a thin accretion disk. Although this paper has not dealt with 
magnetically-driven outflows, these are not incompatible with the HH-driven out flows. The 
hypothesis of massive star formation presented here, following iKetol (j2002bl . 120031 ) . suggests 
that both occur. In this hypothesis, high mass stars grow by accretion from lower mass stars. 
Since magnetically-driven outflows are an inescapable part of low-mass star formation, the 
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massive stars-to-be must reach spectral type B with a magnetically-driven outflow. 



As the star continues to grow to earlier type B, and the i onizing radiation creates an 
HII region, what becomes of the magnetically driven outflow? iBeuther fc Shepherd! (120051 ) 
suggest that there is a of lack observations of collimated jet-like outflows associated with stars 
earlier than type Bl compared with the number of observations of jet-like outflows from later 
type B stars. At the moment we do not know whether the jets are destroye d by processes 
assoc iated with the ionization or survive the formation of the HII region (ITan fc McKee 



20031 ) but are difficult to detect. 



9. Comparison with observations 



A model of accretion with a bipolar ionized outflow has been compared to observations 
of both ionized gas a nd molecular gas in and around the ultracompact HII region G10.6-0.4 



(IKeto fc Wood! 120061 ). The model for the accretion flow in G10.6-0.4 is based on the same 
model of streamlines on ballistic trajectories presented here. These observations are unique 
in mapping the velocities in both the molecular and ionized gas and in following the accretion 
flow from the molecular phase through the ionized pha se. However, the s ource G10.6-0.4 is 



not unique. Molecul ar line observations of G28.20-0.05 flSollins et al.ll2005f ) and G24.78+0.08 
(IBeltran et al.l 120061 ) have also been interpreted as accretion disks and outflows. Other 
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10. Conclusions 



This paper provides a simple theoretical description of some of the effects of the ioniza- 
tion of a massive-star forming accretion flow. The model presented here is based o n three 
simp le models: (1) accretion with rotation as streamlines o n ballistic tr ajectories (jUlrich 
1976 ); (2) an outflow as a pressure-driven isothermal wind (Parker 1958h ; and (3) an ion- 
ization front gravitationally trapped within an accretion flow ketoll2002bh . The composite 
model applies to the evolutionary phase in the formation of massive stars of type B and 
earlier when the stars are both hot enough to ionize an HII region around the star and yet 
are still growing by accretion. The model accretion flows and HII regions may be described 
in terms of three characteristic radii which completely determine the morphology and flow 
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pattern. Variation of the relative magnitudes of the parameters suggests how different mor- 
phologies may be related to the same underlying model. The models further suggest an 
evolutionary sequence driven by the increasing ionization of a star growing by accretion. 
These hypercompact to ultracompact HII regions are different in that their dynamics are 
dominated by gravitational force of their stars whereas the larger HII regions are dominated 
by the thermal pressure of the ionized gas. 
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Fig. 1. — Model of a high angular momentum accretion flow subject to 3 levels of ionizing 
radiation, low {top), medium (medium), and high (bottom) as defined in §6j The figures show 
the log of the density of molecular gas in blue (top) and of the ionized gas in red (middle 
and bottom) in a slice in the XZ plane of the flow. The color scales range from to 1.6 x 10 7 
cm -3 (molecular) (top), from to 1.2 x 10 7 cm -3 (ionized) (middle), and to 1.3 x 10 7 cm -3 
(ionized) (bottom). The circle shows the location of the Bondi-Parker critical radius of the 
ionized gas for spherical flow. The arrows show the velocity of the flow in the XZ plane. In 
the top figure, the longest arrow in the molecular flow represents 26.6 kms -1 and the longest 
arrow in the ionized flow represents 21.5 kms -1 . In the middle figure, the longest arrow in the 
molecular flow represents 8.0 kms -1 and the longest arrow in the ionized flow represents 28.2 
kms -1 . In the bottom figure, the longest arrow in the molecular flow represents 5.4 kms -1 
and the longest arrow in the ionized flow represents 29.4 kms -1 . In the ionized outflow flow, 
the velocity is the sound speed at the critical radius. The axes are labeled in units of Rd, 42 
AU (left) and 34 AU (right). 



-18- 




Fig. 2. — Two accretion flows that differ in the relative amount of angular momentum. The 
flow in the right column has 2 times the angular momentum of the flow in the left column. 
The two figures show a slice in the XZ plane with the log of the density of the molecular 
gas in color and the velocity as arrows. The density ranges from 4.5 x 10 3 to 4.5 x 10 6 cm -3 
(molecular) (left) and from 4.5 x 10 3 to 5.2 x 10 6 cm -3 (right). The longest arrow in the 
molecular flow (blue arrows) represents 1.8 kms -1 (left) and 0.4 kms -1 (right). The yellow 
arrows show the velocity of the ionized outflow at half the scale. 



